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Abstract 

Investigating the long time asymptotics of the totally asymmetric 
simple exclusion process, Sasamoto obtains rather indirectly a formula 
for the GOE Tracy- Widom distribution. We establish that his novel 
formula indeed agrees with more standard expressions. 



1 Introduction 

The Gaussian orthogonal ensemble (GOE) of random matrices is a probabil- 
ity distribution on the set of N x N real symmetric matrices defined through 

Z is the normalization constant and dH = rii<j<j<7v d-ffij. The induced 
statistics of eigenvalues can be studied through the method of Pfaffians. Of 
particular interest for us is the statistics of the largest eigenvalue, E\. As 
proved by Tracy and Widom [H], the limit 

lim P(E l <2N + sN 1/3 ) = F^s) (2) 

exists, P being our generic symbol for probability of the event in parenthesis. 
Fi is called the GOE Tracy- Widom distribution function. Following [3| it can 
be expressed in terms of a Fredholm determinant in the Hilbert space L 2 (H) 
as follows, 

F^s) 2 = det (1 - P S (K + \g)(f\)P s ), (3) 
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where K is the Airy kernel defined through 

K(x,y) = [ dX Ai(x + A) Ai(y + A), 

g(x) = Ai(x), (4) 
f(y) = 1- [ dX Ai(y + A), 

and P s is the projection onto the interval [s, oo). 

The GOE Tracy- Widom distribution Fi(s) turns up also in the theory of 
one-dimensional growth process in the KPZ universality class, KPZ standing 
for Kardar-Parisi-Zhang jlj. Let us denote the height profile of the growth 
process at time t by h(x, t), either x G R or x G Z. One then starts the growth 
process with flat initial conditions, meaning h(x, 0) = 0, and considers the 
height above the origin x = at growth time t. For large t it is expected 
that 

h(0,t) = Cl t + c 2 t 1/3 £i. (5) 

Here C\ and c 2 are constants depending on the details of the model and £i is 
a random amplitude with 

P(£i<s)=F 1 (8). (6) 

For the polynuclear growth (PNG) model the height h(0, t) is related 
to the length of the longest increasing subsequence of symmetrized random 
permutations for which Baik and Rains [Tj indeed prove the asymp- 
totics ©, @, see for further developments along this line. Very recently 
Sasamoto P succeeds in proving the corresponding result for the totally 
asymmetric simple exclusion process (TASEP). If rjj(t) denotes the occupa- 
tion variable at j G Z at time t, then the TASEP height is given by 

f2iV t + Eii(l-2^(t)) forj>l, 
h(j,t) = hN t forj=0, (7) 

l2tft-Xt i+ i(l-2»*(t)) forj<-l, 

with N t denoting the number of particles which passed through the 
bond (0, 1) up to time t. The flat initial condition for the TASEP 
is ...010101.... For technical reasons Sasamoto takes instead 
...010100000... and studies the asymptotics of h(—3t/2,t) for large t 
with the result 

h(-3t/2,t) = lt + ^^ SA . (8) 
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The distribution function of the random amplitude £sa is 



P(£sa <s) = F SA (s) 



(9) 



with 




determinant is in L 2 (R). 

The universality hypothesis for one-dimensional growth processes claims 
that in the scaling limit, up to model-dependent coefficients, the asymptotic 
distributions are identical. In particular, since © is proved for PNG, the 
TASEP with flat initial conditions should have the same limit distribution 
function, to say 



As written above, the s-dependence sits in the projection P s . It will turn out 
to be more convenient to transfer the s-dependence into the integral kernel. 
From now on the determinants are understood as Fredholm determinants in 
L 2 (R + ) with scalar product (-, ■). Thus, whenever we write an integral kernel 
like A(x, y), the arguments are understood as x > and y > 0. 
Let us define the operator B(s) with kernel 



By |Zj ||-B(s) 2 || < 1 and clearly B(s) is symmetric. Thus also ||S(s)|| < 1 
for all s. B(s) is trace class with both positive and negative eigenvalues. 
Shifting the arguments in ()1U|) by s, one notes that 




(11) 



2 The identity 



B(s)(x, y) = Ai(x + y + s). 



(12) 



F SA (s) = det(l-B(s)). 



(13) 



Applying the same operation to © yields 

Fx(s) 2 = det (1 - B{sf 



9)(f\) 



(14) 



with 



M 




Ai(x + s) = (B(s)5)(x), 
1- / d\Ai(y + X + s) 




(15) 
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Here 5 is the 5-function at x = and 1 denotes the function l(x) = 1 for all 
x > 0. 5 and 1 are not in L 2 (R + ). Since the kernel of B(s) is continuous 
and has super-exponential decay, the action of B(s) is unambiguous. 

Proposition 1. With the above definitions we have 

det(l-B(s))=Fi(s). (16) 

Proof. For simplicity we suppress the explicit s-dependence of B. We rewrite 

F^sf = det((l-B)(l + B-\B5)(l\)) 

= det(l - B) det(l + B) (l - (5, B(l + 

= det(l- B)det(l + B)(5,(l + ByH) (17) 

since 1 = (5, 1). Thus we have to prove that 

det(l - B) = det(l + B)(6, (1 + B)' 1 !}. (18) 

Taking the logarithm on both sides, 

lndet(l - B) = lndet(l + B) + In (5, (1 + (19) 

and differentiating it with respect to s results in 

- 1* ((1 - B)-f B B)) = IV ((1 + B)-!*)) + f^^f (20) 
where we used 

A M det(T)) = Tr(T- 1 ^T). (21) 

Since B(s) — > as s — > oo, the integration constant for (f2T)|) vanishes and we 
have to establish that 



Define the operator D = Then using the cyclicity of the trace and 
Lemma 

-2Tr((l-£ 2 ) -1 — £)) = -2Tr ((1 - B^DB)) 

= (5,(1- B 2 Y X BS). (23) 
Using Lemma El and Dl = 0, one obtains 

(6, j- g (t + B)-H) = (6, (1 - B 2 y 1 B5)(S, (1 + BY 1 !). (24) 
Thus d21 follows from 1)25 ) and l|25 J) . □ 
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Lemma 2. Let A be a symmetric, trace class operator with smooth kernel 
and let D = j- . Then 

2 Tt(DA) = —(5, A5) (25) 
where DA is the operator with kernel -^-A(x,y). 

Proof. The claim follows from spectral representation of A and the identity 
dxf'(x)f(x) = -/(0)/(0) - / dxf(x)f(x). (26) 

□ 

Lemma 3. It holds 

J?-(i + 5)- 1 = (1 - B 2 Y X BD + (1 - P 2 )- 1 ^)^! + P)" 1 !. (27) 
as 

Proo/. First notice that §-B = B = DB. For any test function /, 
(Bf)(x) = [ dyd y Ai(x + y + s)f(y) 

= -Ai(x + s)/(0)- / dyAi^ + y + s)/'^). (28) 

Thus, using the notation P = \B5)(5\, one has 

DP = -BD - P. (29) 
Since ||P|| < 1, we can expand + P) _1 in a power series and get 

r) Pi 71—1 

^-(1 + P)- 1 = y^(-i)"^-S» = V(-l)" V B k DB n ~ k . (30) 
as as / — ' 

n>l n>l k=0 



Using recursively ()29|) we obtain 



n— 1 / \ n n— 1 n— 1 

^ B k DB n-k = "I" j j g»^ + ^^(_l)i+l j gfcp^"-fc-l 

fc=0 j=0 k=j 

n-1 „ , / -, \£ 



' ! ""p"P + ^ 1 + ( 2 X) B k PB n - k -\ (31) 



fc=0 

Inserting (|31jl into (jBTlj) and exchanging the sums results in 

^(1 + P)" 1 = ^P 2n+1 P + ^ ^ 1 + ( ~ 1)fc p fc P(-P)"-( fc+1 > 

n>l fc>0 n>fc+l 

= (1-P 2 )" 1 PP+(1-P 2 )" 1 P(1 + P)" 1 . (32) 

□ 



3 Outlook 



The asymptotic distribution of the largest eigenvalue is also known for Gaus- 
sian unitary ensemble of Hermitian matrices (/3 = 2) and Gaussian symplectic 
ensemble of quaternionic symmetric matrices (/3 = 4). As just established, 
for = 1, 

Fx(s) = det(l - B(s)), (33) 

and, for f3 = 2, 

F 2 (s) = det(l- J B(s) 2 ), (34) 

which might indicate that F^s) equals det(l — B(s) 4 ). This is however 
incorrect, since the decay of det(l — _B(s) 4 ) for large s is too rapid. Rather 
one has 

F A (s/V2) = -(det(l-5(s)) + det(l + £(»)). (35) 

This last identity is obtained as follows. Let U(s) = | f s °° q(x)ds with q the 
unique solution of the Painleve II equation q" = sq + 2q 3 with q(s) ~ Ai(s) 
as s — > oo. Then the Tracy- Widom distributions for {3 = 1 and (3 = 4 are 
given by 

Fx(s) = exp(-f/(s))F 2 (s) 1 / 2 , F 4 (s/v^) = cosh(f/(s))F 2 (s) 1 / 2 , (36) 
see jEj. Thus F 4 (s/V2) = |(Fi(s)+F 2 (s)/F 1 (s)), from which ® is deduced. 
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